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Abstract 

The influence of the Coulomb distortion for quasi-elastic (e,e') scattering on highly 
charged nuclei is investigated in distorted wave Born approximation for electrons. 
The Dirac equation is solved numerically in order to obtain exact electron continuum 
states in the electrostatic field generated by the charge distribution of an atomic 
nucleus. Different approximate models are used to describe the nucleon current in 
order to show that, at high electron energies and energy- momentum transfer, the 
influence of Coulomb distortions on (e,e') cross sections can be reliably described 
by the effective momentum approximation, irrespective of details concerning the 
description of the nuclear current. 
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1 Introduction 



Quasi-elastic electron scattering on nuclei has represented during nearly three 
decades one of the most successful tools to study the nuclear and nucleon 
structure. Experiments have been performed at SLAC [1,2,3], the MIT Bates 
Laboratory [4,5,6,7,8,9,10,11,12], the Saclay Laboratory [13,14,15,16,17,18], 
and more recently at JLab in order to explore this reaction. Both inclusive 
(e, e') and exclusive (single arm (e, e'N) or double arm (e, e'NN)) experiments 
contributed to a deeper understanding of the many-body structure of strongly 
interacting systems like light and heavier nuclei opening the possibility of 
investigating also the in-medium nucleon properties. 
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Inclusive (e, e') scattering, where only the scattered electron is observed, pro- 
vides information on a number of interesting nuclear properties like, e.g., the 
nuclear Fermi momentum [19], high-momentum components in nuclear wave 
functions [20], modifications of nucleon form factors in the nuclear medium 
[21,22], the scaling properties of the quasi-elastic response allow to study the 
reaction mechanism [23], and extrapolation of the quasi-elastic response to 
infinite nucleon number A = oo provides us with a very valuable observable 
of infinite nuclear matter [24]. 

However, experimental studies of inclusive and exclusive reactions induced 
by electrons are hampered in the case of target nuclei with a large number 
of protons by the strong Coulomb field, which induces a distortion of the 
electron wave front, hereby modifying the structure of the (e, e') cross section 
and inducing sizable effects in the longitudinal and transverse separation of 
the electromagnetic response [25,26,27,28,29]. 

Still, there is a considerable theoretical and experimental interest in extracting 
longitudinal and transverse structure functions as a function of energy loss for 
fixed three momentum transfer for a range of nuclei. One of the important 
topics mentioned above concerning the in-medium modification of the nucleon 
form factors is related to the question whether the Coulomb sume rule is 
violated in nuclei. This sum rule states that the number of protons in a nucleus 
can be obtained from an integral of the electric response of the nucleus over the 
full range of the electron energy loss at large three-momentum transfer [22]. 
Unfortunately, the conclusions reached by the large experimental programs 
carried out at Bates, Saclay and SLAC for a variety of nuclei ranged from a 
full saturation of the Coulomb sum rule to its violation by 30%. 

The approval of the Thomas Jefferson National Accelerator Facility (TJNAF) 
Proposal E01-016, entitled 'Precision measurement of longitudinal and trans- 
verse response functions of quasi-elastic electron scattering in the momentum 
transfer range 0.55 GeV < \q\ < 1.0 GeV, resulted in experiments performed 
very recently in Hall A at TJNAF, using target nuclei like 4 He, 12 C, 56 Fe, 
and 208 Pb. In the present study, we therefore focus on the particular case of 
208 Pb(e,e') scattering. Preliminary calculations for 56 Fe(e, e') show that the 
conclusions drawn below in the present study apply to heavy nuclei like 56 Fe 
in an analogous manner as well. 

Since the plane wave Born approximation (PWBA) for electrons is no longer 
adequate for the calculation of scattering cross sections in the strong and long- 
range electrostatic field of highly charged nuclei, it has become clear in recent 
years that the correct treatment of the Coulomb distortion of the electron 
wave function due to the electrostatic field of the nucleus is unavoidable if 
one aims at a consistent interpretation of experimental data. The theoretical 
framework to investigate Coulomb corrections to the electron-nucleus cross 
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sections is well established and is called distorted wave Born approximation 
(DWBA) in contrast to the better known PWBA where the incoming and 
outgoing charged leptons are described by (Dirac) plane waves, neglecting the 
effect of the Coulomb interaction between the projectile and the target nucleus. 
The application of the DWBA scheme is in principle straightforward, but the 
numerical complications lead to extremely time consuming calculations. 

DWBA calculations with exact Dirac wave functions have been performed 
by Kim et al. [30] in the Ohio group and Udias et al. [31,32] for inclusive 
quasi-elastic scattering on heavy nuclei. However, these cumbersome calcula- 
tions are difficult to control by people who were not directly involved in the 
development of the respective programs. Early DWBA calculations for 12 C 
and 40 Ca were presented in [26]. Coulomb corrections have also been evalu- 
ated theoretically by a group from Trento University [33], where it was found 
that the standard method (used mainly in the case of light nuclei) to han- 
dle Coulomb distortions for elastic scattering in data analysis, the so-called 
effective momentum approximation (EMA), works with an accuracy better 
than 1% for the description of Coulomb distortion effects, if the energy of the 
scattered electrons and the momentum transfer is sufficiently high. On the 
other hand, the Ohio group derived significant corrections beyond the EMA. 
The findings presented in this paper confirm that the conclusions drawn by 
the Trento group are basically correct. There is a certain mismatch between 
the EMA and exact calculations for small energy transfer, which, however, 
strongly decreases when the energy transfer becomes significantly larger than 
the typical removal energy of the bound nucleons. 

Calculations using an eikonal approximation (called eikonal distorted wave 
Born approximation, EDWBA) for electrons have been performed by the Basel 
group [34,35], which seemed to confirm the results of the Ohio group. However, 
a poor approximation for the focusing of the electrons near the nucleus was 
used, leading to an overestimation of the Coulomb corrected cross sections. 
Through a subsequent analysis using exact electron wave functions we found 
that the eikonal approximation with an improved non-perturbative descrip- 
tion of the electron wave function amplitude supports the observation that the 
EMA is indeed a valid tool for the description of Coulomb distortion effects 
[36]. A similar preliminary result was obtained in [37], where exact solutions 
of the Dirac equation were used for the electron wave functions. The nuclear 
current was modeled within the framework of a simplified single particle shell 
(SPS) model, with harmonic oscillator wave functions for the bound nucleons 
and plane waves for the knocked-out nucleons. The cross section for inclu- 
sive quasi-elastic scattering process is usually calculated approximately by 
summing over all knock-out processes, where the individual bound nucleons 
inside the nucleus leave the nucleus after a sufficiently high energy transfer 
uj = €i — €f from the scattered electron, where q is the initial and e/ the fi- 
nal electron energy. The recoil nucleons move in an energy dependent optical 
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potential. It is common practice that the imaginary part of the optical po- 
tential is not taken into account in the calculations, since the imaginary part 
is intended to describe the loss of nucleon flux inside the nuclear medium. In 
inclusive processes, only the electrons are observed, and one may argue that 
it is not necessary to take into account whether or not the recoil nucleons 
actually leave the nucleus or initiate some subsequent nuclear reactions. Still, 
it has been shown in [38] that (e, e') peak intensities are reduced by typically 
5 — 10% as a result of final state interaction (FSI) for the kinematics relevant 
for this paper, and how this reduction is related to the imaginary part of the 
final state nucleon optical potential (see eq. (29a) in [38]). The effects of the 
FSI at high momentum transfer have been discussed in [39]. 

In this paper, we present calculations based on two different models for the 
nuclear current. The final state nucleons are not described by wave functions 
obtained as solutions of the Dirac equation with some nuclear model poten- 
tial. Instead, an eikonal approximation and the plane wave approximation are 
used, reducing the computational costs in a very effective way. It turns out 
that the use of approximate wave functions is not a weakness, since we find 
that at sufficiently high energies and momentum transfer, the effective momen- 
tum approximation provides an accurate description of Coulomb corrections, 
irrespective of the nuclear model used. It is a fact that also elaborate nuclear 
current models, based, e.g., on solutions of the Dirac equation for the nucleons 
in some relativistic nuclear model, cannot be considered as a fully satisfying 
strategy, as the still sizeable discrepancy between measured and calculated 
cross sections shows (see, e.g., [30]). Computational costs always restrain us 
from using more realistic nuclear models, where effects like nucleon correla- 
tions and meson exchange currents can be properly taken into account. The 
electrons are described by exact solutions of the Dirac equation in our DWBA 
calculations. 



2 From DWBA to EMA 



Inclusive (e, e') cross sections are usually calculated by integrating over all 
final state nucleons in the knock-out (e, e'N) cross sections obtained from a 
single-particle shell model [40] for the bound protons and neutrons, as will 
be explained in further detail in the following section. The differential cross 
section for single nucleon knockout is given by [31] 



dejdVtfdEfdVtf 
Aa 2 



e 2 f E fPf S(e l + E A -e f -E f -E A ^)Y: \W tf \ 2 , (1) 
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with the matrix element 

'{r e -r N ) 



W if = j^y 3 jd 3 r e j #r N ]W\%<fe) q/2+l0 (2) 



where j%{f^) is the nucleoli transition current obtained within the framework 

of some suitable nuclear model, the J2 in eq. (1) indicates the sum (average) 

over final (initial) polarizations, and Ea, Ea-i is the energy of the initial and 
final nucleus, respectively. The virtuality of the exchanged photon is q' 2 = 
uj 2 — q' 2 , where q° = uj = — e/ is the energy tranferred by the electron to 
the nuclear system. In our calculations, the final nucleus was considered as 
a spectator, and the final energy Ef of the nucleon was calculated from the 
energy transfer reduced by a reasonable separation (or removal) energy. 

In the PWBA, the electron current is given by 

#(r.) = u^TX^e**-*'*, (3) 



where u Si ,u Sf are initial/final state plane wave electron spinors corresponding 

to the initial/final electron momentum kij and spin Sij, normalized according 
to 



u s (k)^°u s >(k) = 5 SyS >, 



(4) 



and a short calculation yields a three-dimensional integral 



(5) 



with the four-momentum transfer squared Q 2 = (fcj — kf) 2 — oj 2 = q 2 — u 2 . 

In the DWBA, exact solutions of the Dirac equation for electrons in the elec- 
trostatic field of the nucleus are used instead of plane wave functions. In a 
general sense, the EMA accounts for the two effects of the Coulomb distor- 
tion, namely the momentum enhancement of the electron near the attractive 
nucleus and the focusing of the electron wave function. Therefore, transition 
amplitudes are calculated by replacing the asymptotic initial and final state 
momenta fcj, kf of the electron by appropriate effective momenta in the PWBA 
transition amplitude. Additionally, the cross section obtained from the effec- 
tive transition amplitude is multiplied by an appropriate focusing factor, which 
accounts for the fact that the electron wave function amplitude is enhanced in 
the vicinity of the nucleus, i.e. from a semiclassical point of view, the nucleus 
acts like a lens and concentrates the electrons towards the nuclear interior. 
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Indeed, one may establish a direct connection between the DWBA and the 
EMA. The DWBA transition amplitude comprises the exact solution of the 
Dirac equation for the electron in the electrostatic field of the nucleus. In this 
way, the DWBA accounts for the exchange of the 'soft' photons between elec- 
tron and nucleus, whereas the explicit photon propagator in eq. (2) accounts 
for the exchange of one hard photon. 

Using the distributional identity 

r d 3 q' _ e ±iu}r 

J (27r)3 w2 _ g -,2 ± . --4^T> U 



one obtains 



W if = -^J d 3 r e d 3 r N j»((r e )tf(r N ) 



r e - r N 



1 r (? ujr 
47T J d3red3rN {p e ^P N ^ N ^ ~~ je(r e )jN(r N )} — , 

r = \f e - f N \, (7) 



with p e , p^, j e , jw denoting the transition charge and current densities of the 
electron and the nucleon, respectively. The double volume integral presents a 
clear numerical disadvantage of this expression. According to Knoll [41], one 
may introduce the scalar operator 




such that the transition amplitude can be expanded in a more convenient form 
1 



W tf = 



d A r 



p N (r)Sp e (r) -j N Sj e (r) 



(9) 



The single integral is limited to the region of the nucleus, where the nuclear 
current is relevant. The expansion eq. (8) is an asymptotic one, in the sense 
that there is an optimum number (depending on q) of terms that give the 
best approximation to the exact value. Considering terms up to second order 
in the derivatives only one obtains 



i + 



2iq V + A 4(<f V)' 



-iqr 



Q 2 



(Q 2 ¥ 



Pe{r) 
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+ current terms 



(10) 



This approximation has been used in our calculations [37] and in an equivalent 
way in [30] , where it was also observed that higher order terms can indeed be 
neglected for the calculations relevant for this work. It should also be men- 
tioned that the approximation has been applied for the first time to inclusive 
reactions in [27,28], and discussed in the context of the (e,e'p) quasi-elastic 
reaction in [42]. 

There is indeed a simple semiclassical interpretation of this expansion. If 
one assumes that the electron is scattered locally with a momentum trans- 
fer qiocif) = Q + Aq(r), calculable from the classical local momentum of the 
electron in the electrostatic potential and which deviates by Aq from the the 
asymptotic momentum transfer q, then one has 

Q 2 _ q 2 ~UJ 2 Q 2 
QL (<? + MY Q 2 + 2qAq + Aq 2 " 

2qAq + Aq 2 4(qAq) 2 

Q 2 (Q 2 ) 2 1 } 



Furthermore, one may note that the PWBA electron transition current can 
be written as 



Applying the operator q = j on j£(r) leads to qj£(r) = qje(r), i.e. j acts as 
momentum transfer operator and replacing q in (11) by q, one immediately 
obtains the differential operator part of S. 

Therefore, the Knoll operator S first factors out the undistorted e iqr -term 
in the current, then calculates the ratio of the asymptotic momentum trans- 
fer with a quantity that can be interpreted as the local momentum transfer 
generated from the distorted electron current, and finally reinserts the e~ iqr - 
behavior of the undistorted current. 

Assuming that one can approximate the effect of the S'-operator by replacing 
it by an average (effective) momentum transfer Qlff, one obtains as a first 
approximation for the transition matrix element 

W if = / JrjStftftf, (13) 
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keeping in mind that the DWBA electron current still has to be used in eq. 
(13). Since a highly relativistic particle (with negligible mass m) is moving 
nearly on a straight line, the classical momentum of a particle with asymptotic 



momentum k and energy y k 2 + m? ~ \k\ moving inside a potential V(r) is 
given by -^(k — V(r)). Additionally, since the knockout process is nearly local 

for large momentum transfer (e.g., Q 2 = (400 Me V) 2 corresponds to a photon 
propagation length scale of ~ \ fm which is much smaller than the extension of 
a heavy nucleus), and since it takes place in the entire volume of the nucleus, 
it is self-evident that one should calculate Qlff from a potential value V which 
is obtained by an averaging process over the nuclear density profile p(r) 

= (k? f - kff - c 2 , k$ = ^L(k hf - V) (14) 



i,f\ 



with 



V = l^P^, (15) 
J d 6 rp(r) 



where p can also be used for an approximate description of the charge density 
of the nucleus. In the case of a homogeneously charged sphere with radius 
R and charge number Z, the electric potential in the center of the sphere is 
given by V(0) = — |^ (a is the fine structure constant), whereas the potential 
averaged over the volume of the sphere is given by V — 4V r (0)/5. 

A further approximation which can be made in eq. (13) stems from the ob- 
servation that the attractive nucleus focuses the electron wave function in the 
nuclear region. It can be shown that for highly relativistic energies, the am- 
plitude of the initial/final state electron wave function in the nuclear center 

— * 

is enhanced by a focusing factor fij(0) = (1 — V(0)/\kij\), however, solving 
the Dirac equation for an electron in the electrostatic field of a highly charged 
nucleus again shows that the focusing averaged over the nuclear volume is 
again well described by fij = (1 — V/\kij\). Since the electron transition cur- 
rent is given by j£ = vpej^tpe^, one may replace j£(r) by fi^ ff^j£' EMA , 
where j£ EMA is calculated with plane electron waves with effective momenta. 
A detailed discussion of the focusing effect can be found in [43] . 

If the same average potential value V is used to calculate the effective momenta 
and the effective focusing factors, then a short calculation shows that the 
focusing factors cancel against the effective photon propagator ~ Q7ff ^ n ^ ne 
cross section in the sense that 



Uf f )VfY i 
Qltf Q 



j ■ V (16) 
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and consequently, the EMA cross section is calculated by replacing the electron 
current in (10) by the PWBA current calculated from the effective momenta, 
and leaving everything else unchanged. A different, but completely equivalent 
strategy is to calculate first the theoretical cross section based on the effec- 
tive momenta and energies instead of the corresponding asymptotic values. 
The intermediate result obtained this way must be multiplied subsequently 
by (fi^) 2 in order to account for the focusing of the incident electron. The 
focusing of the final state electron is automatically taken into account by the 
artificially enhanced phase space factor {e e ^) 2 ~ \k e ^\ 2 . 

The cancellation mechanism that plays between the focusing effect and the the 
photon propagator in the matrix element is accidental and not really exact. 
It works in a very satisfactory way for nuclear charge distributions which 
are close to the homogeneous, spherical case. If the energy of the final state 
electron becomes significantly smaller than 150 MeV and the four-momentum 
transfer is below (300 MeV) 2 , the semiclassical picture on which the EMA is 
based starts to fail. Furthermore, for small energy transfer oo comparable to 
the removal energy of the nucleons, details of the nuclear structure and the 
interaction of the recoil nucleons with the nuclear matter become increasingly 
important, which makes the applicability of the EMA questionable. However, 
this kinematical region is of minor interest in future experiments. 



3 Models and approximations 

3. 1 Single particle shell model 

As mentioned before, inclusive (e, e') cross sections are commonly calculated 
by integrating over all final state nucleons in (e, e'N) cross sections obtained 
from a single-particle shell (SPS) model for the bound protons and neutrons. 
Usually, the fact that there are short-range and tensor correlations between 
nucleon pairs which lead to a partial depletion of the single particle shells is 
ignored in SPS (e, e') calculations. Within a correlated basis functions theory 
calculation, it was found in [44] for nuclear matter that there is an approximate 
probability of 20% for a nucleon being in an correlated state [45]. However, 
correlated nucleons account for approximately 37% of the average removal en- 
ergy. In view of the fact that the average binding energy of, e.g., protons is 
21 MeV in our SPS calculations, this observation implies that the correlated 
nucleons correspond to strongly bound particles with binding energies around 
100 MeV, if an average removal energy of 37 MeV is assumed. The impact 
of correlated nucleons on the inclusive cross section is therefore strongly sup- 
pressed for an energy transfer clearly smaller than 150 MeV compared to the 
strength of the uncorrelated nucleons. Still, a reduction of ~ 20% of the peak 
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strength results from the suppression of correlated nucleons, combined with 
a slight shift of the peak towards higher energy transfer 00 and an increase of 
strength in the high-u; tail. 



A further observation, which is of minor importance concerning the present 
calculations, is the fact that the rms radius of uncorrelated nucleons tends to 
be slightly larger than the total rms radius of the nucleus [45]. We therefore 
performed calculations within a strongly simplified framework, where the cor- 
related (high-momentum) nucleons where described by model wave functions 
with a high-momentum component and large binding energies. We found that 
even in this case, the EMA remains valid to an acceptable extent, if the 'cor- 
related' nucleons are distributed in a reasonable homogeneous manner inside 
the nucleus, and if the momentum transfer and final electron energy are suf- 
ficiently large. In this work, we therefore focus on the main contribution of 
the uncorrelated nucleons in the region of the quasi-elastic peak, in order to 
facilitate a comparison to other works [30]. 

The relativistic form of the wave functions for bound protons and neutrons was 
constructed from the non-relativistic wave function. Non-relativistic bound nu- 
cleon wave functions were generated using a self-consistent Schrodinger solver 
for a non-relativistic Woods-Saxon potential and the additional Coulomb po- 
tential for the protons. The Woods-Saxon potential, which included an LS 
coupling term, was optimized in such a way that the experimental binding 
energies of the upper proton shells and the rms of the nuclear charge dis- 
tribution were reproduced correctly. The standard way to construct nucleon 
four-spinors ipN from a Schrodinger-type nucleon wave function xn was then 
applied, given by 



where p denotes the differential momentum operator, E is the energy of the 
nucleon, and N is a normalization factor which is typically close to one. The 
same strategy, which is exact in the case of Dirac plane waves, was also applied 
to construct the small component of the Dirac spinor in the case of the nucleon 
eikonal approximation. For a discussion of this approximation we refer to [46]. 

The nuclear charge distribution, which is the main input for the calculation 
of the Dirac electron wave functions, and nucleon optical potentials of heavy 
nuclei are often approximated by the help of a Woods-Saxon distribution. Since 
the normalization and moments of Woods-Saxon distributions are rarely found 
in the literature, some important expressions are listed in Appendix A. 




(17) 
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3.2 Eikonal approximation for nucleons 



In order to generate realistic nucleon wave functions in a computationally 
efficient way, we used eikonal wave functions for the final state protons and 
neutrons, with the aim to compare the (Coulomb corrected) cross sections 
obtained from the eikonal nuclear model to the cross section obtained from 
simple nucleon plane wave calculations. This comparison will serve as a test 
whether the applicability of the EMA is strongly dependent on the nuclear 
model or not. 

The eikonal phases xi,2 generated by a potential U(f) take the local (semi- 
classical) momentum modification of the incoming/outgoing nucleon with ini- 
tial/final asymptotic momentum pij = \Pij\Pij approximately into account 
by modifying the plane waves describing the nucleons according to 

e Wi,ff _^ e mjr±ixi,2(r) 



where 



Xi(r) = -prr / U(f+ Pis)ds , z = far ', (19) 

\Pi I J 



or 



Xi(r) = -— \ f U(x,y,z')dz' , z = far (20) 



if the initial nucleon momentum is parallel to the z-axis. The eikonal phase 
for the final state relevant for the calculations presented in this work is 

F °r 

X2{?) = -^^U{r+p f s')ds' , (21) 

where E i: f = (pfj + m 2 N ) l l 2 is the energy of the nucleon with mass m N , 
and Ef/\pf \ is the velocity of the knocked-out nucleon. An energy-dependent 
volume-central part of an optical Woods-Saxon type model potential V nuc (r, E) 
as given in a recent work [47] was used for our calculations, as well as the 
energy-independent Woods-Saxon potential which was used to generate the 
bound state nucleon wave functions. 

The total potential U(r,E) = V nuc (r,E) + V cou i(r) is given as the sum of a 
Woods-Saxon potential and the Coulomb potential of the final state nucleus 
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in the case of protons. The depth of the Woods-Saxon type nuclear potential 
V nuc (r,E) = - - V Y S{E 1 (22) 



depending on the energy E of the proton (all quantities in MeV) is given by 
V WS (E) = Vl [l- v 2 (E - E f ) +v 3 (E- E f f - v 4 (E - E f f] , (23) 

where Ef = —5.9 and 

v l = 67.2, v 2 = 7.9 • 10" 3 , v 3 = 2.0 • 10" 5 , v 4 = 7 • 10" 9 , (24) 

whereas the parameters for neutrons are Ef = —5.65 and 

Vl = 50.6, v 2 = 6.9 • 10" 3 , v 3 = 1.5 • 10" 5 , v A = 7 ■ 10~ 9 . (25) 

The range and diffusivity are both for protons and neutrons given by r nuc = 
7.371 fm and a nuc = 0.646 fm according to [47]. 

The imaginary part of the optical potential, which is intended to describe the 
loss of flux in proton (neutron) elastic scattering, has been neglected in our 
(e, e') calculations. As mentioned in the introduction, in inclusive processes, 
only the electrons are observed, and the eventuality whether or not an ejected 
nucleons got 'lost', i.e. initiated some subsequent nuclear reaction, does not 
play a very important role. Furthermore, the FSI affects the inclusive cross 
sections in a similar way whether the electron wave function distortion is taken 
into account or not, such that a comparison of electron PWBA and DWBA 
calculations is still possible. 

To calculate the eikonal integral for a Woods-Saxon potential numerically 
would be relatively time consuming, and unfortunately, the eikonal integral 
of a Woods-Saxon potential is not suitable for simple analytic estimations 
[48,49]. Therefore, one may approximate the Woods-Saxon shape of the mean- 
field optical nucleon potential V nuc in an effective way by a power series in the 
distance from the nuclear center. We found that the following simple two- 
parameter approximation provides an efficient approximation for V nuc (r) in 
the range < r < R max in the case of 208 Pb, where R max = 10 fm: 

V nuc (r, E) = -V WS {E) + Cl (E)r 7 + c 2 (E)r 6 (26) 



with 



Ci = - 6 bl~j 02 = 7 wr- (27) 



max max 
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Distance from nuclear center [tm] 



Fig. 1. Comparison of the Woods-Saxon type nuclear potential and the simple poly- 
nomial approximation described in the text, which allows to derive a simple analytic 
expression for the nuclear eikonal phase. 

The approximate nuclear potential V nuc (r, E) satisfies the boundary condi- 
tions V n uc(Rmax,E) = and j;Vnuc(Rmax, E) = 0, which are acceptable for a 
distance of 10 fm from the nuclear center. A comparison of the exact Woods- 
Saxon distribution and approximation eq. (26) is shown in Fig. 1. The contri- 
bution of V nuc (r, E) to the eikonal phase of the nucleons can be derived from 
the integrals 



/ (z 2 + b 2 fdz = \z 7 + h 2 z 5 + b 4 z 3 + b 6 z , 
I Vz 2 + b 27 dz -- 



93 , 6 163, 4 3 25, 2 5 1 7 

b 6 z H H b 2 z 5 + -z 7 

128 192 48 8 



35 



x Vz 2 + b 2 + —b s \og(z + Vz 2 + b 2 ). 
128 



(28) 



In order to be compatible with [47], we used the electrostatic Coulomb po- 
tential generated by an homogeneously charged sphere (Z' = Z — 1, R c = 
1.22CL4 1 / 3 = 7.228 fm) 

V coul (r) = \ Rc } 2 c ' ~ (29) 

: r > R c 



for protons. The eikonal phase can easily be constructed from the expressions 
above, however one should note that the long-range part of the Coulomb 
potential must be regularized and one has (r' 2 = b 2 + z' 2 , r 2 = b 2 + z 2 , z = pif) 

X rVl = -<*Z'§- hm [ / (I - 1 W - ]og(R/Ro) = 

— oo 
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fm _20 -20 fm 

Fig. 2. Typical eikonal phase for outgoing neutrons in a plane through the nuclear 

oo 

center (pf points rearwards to the right). The plot shows / V nuc (f + pfs)ds, where 

o 

Vnucif) is an attractive Woods-Saxon potential with a depth given via Vws(E) = 30 
MeV and corresponding parameters for 208 Pb. 

Ei f R {r + z) \ 

" a i og h^l (30) 



or 

,coul 



xT Ul (r) = 




= aZ'—^- \og(kir — kif) + const., (31) 
N 



which is defined up to a constant phase given by the free parameter Rq. Anal- 
ogously, X2° ul which is relevant for our case is given in a space region where, 
e.g., r > Rc and pff > 0, by 



x ^ (f) = -aZ'— J -\og[ 

\Pf\ \ b 2 J 

= aZ'^f log(k f r + k f r) (32) 



with z = pff and b = r 2 — z 2 . 

Figs. 2 and 3 show typical eikonal phases for neutrons (protons) in (a repul- 
sive Coulomb and) a mean field optical potential for 208 Pb with a depth of, 
e.g., 30 MeV. Note that the maximal depth of the surface in Fig. 2 is given 
approximately by the spatial extension of the Woods-Saxon potential multi- 
plied by the potential depth, i.e. ~ 30 MeV x 14 fm ~ 2.13 he. The eikonal 
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fm 



-20 -20 



fm 



Fig. 3. Typical eikonal phase for outgoing protons in a plane through the nuclear 



center. The plot shows the regularized combined phase 1Z J U(f + pfs)ds with 

o 

U(r) = V nuc (r,E) + V cou i(r) for the same nuclear Woods-Saxon potential as in 
Fig. 2 and the repulsive electrostatic potential for 208 Pb with Z' = 81. 

phases depicted in Figs. 2 and 3 must be multiplied, up to a sign, by a energy- 
dependent factor Ef/\pf\ > 1 in order to get the relevant phase which distorts 
the corresponding nucleon Dirac spinor. 

In addition to the eikonal phase correction, we also took the (de-)focusing 
effect of the nucleon wave functions in the energy-dependent optical potential 
into account. We remark that Baker investigated the second-order eikonal 
approximation for potential scattering in the non-relativistic case [50], finding 
thereby an expression for the focusing factor of continuum Schrodinger wave 
functions. For a spherically symmetric potential V(r), one finds for the central 
focusing factor /(0) of the amplitude of a non-relativistic scalar nucleon wave 
function (see eq. (23) in [50]) 



where pf is the asymptotic momentum and Vf the velocity of the parti- 
cle. Roughly speaking, the approximation is valid if the asymptotic kinetic 
energy of the particle is larger than the depth of the disturbing potential 
m 3> E[ in = Ef — m 3> |V(0)|, and the wave length of the particle ~ 2n/pf 
should be significantly smaller than the extension of the potential. For the 
classical particle momentum in the center of the potential p/(0) one has non- 
relativistically 



/(o) ~ l 



V(0) 



(33) 



2p f Vf ' 



p f (0) = V / 2m(4 n - V(0)) 



= sfh^Ul - V(0)/E[J ~ Pf (l - V(0)/2EL), 



(34) 
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such that 



/(O) ~ 1 



V(0) ^ Pf (0) 



(35) 



2 Pf v f \ Pf 



i.e. it is found that the probability density is enhanced by the square root of the 
ratio of the central and asymptotic momenta £>/(0)/p/, contrary to the result 
/(0) =Pf(0)/Pf — (Ef — V(0))/Efmihe highly relativistic case. One may ask 
how the non-relativistic and the highly relativistic regime are connected. A 
classical relativistic analysis of the particle trajectories shows that the central 
focusing is given by the expression 



which interpolates between the non-relativistic and relativistic regime [51]. 
We additionally solved the Dirac equation for Dirac nucleons in the energy- 
dependent volume-central part of the optical potential used in this work for 
different energies, such that the amplitude of the eikonal phase corrected nu- 
cleon spinors could also be corrected by the corresponding energy- dependent 
amplitude modification. 

The focusing and eikonal phase correction has an impact on the outgoing 
nucleon flux, and both the cross sections calculated in our approach and the 
results found in [30] agree in a different, but acceptable manner with the 
experimental data. An overall phenomenological Perrey factor [52] (correcting 
for the violation of unitarity in our simplified approach) which applies in the 
same manner to the nucleon wave functions irrespective of the electron wave 
functions used, only leads to small corrections to the cross sections of the order 
of a few percent and has been neglected. However, even though the distorted 
nucleon wave function approach presented in [30] is more ambitious than the 
present one, it is unable to reproduce the experimental data on a much higher 
level of precision, since it is still based on a SPS model, which fails to account 
for many physical aspects of the inclusive scattering process. The relevant 
observation presented in this paper is the fact that Coulomb corrections turn 
out to be compatible with the EMA, with no strong dependence on the nuclear 
model used. 

It is clear that the same wave functions for the outgoing nucleon have to be 
used for the electron PWBA and exact DWBA calculations, since we are only 
interested in the role of the electronic part of the Coulomb corrections. As 
mentioned above, for bound nucleons wave functions generated from a non- 
relativistic Woods-Saxon potential and a Coulomb potential for the protons 
were used. The Woods-Saxon potential, which included an LS coupling term, 




(36) 
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was optimized in such a way that the experimental binding energies of the 
upper proton shells and the rms of the nuclear charge distribution were repro- 
duced correctly. The rms radius of the neutron density distribution was taken 
from a recent study presented in [53], the binding energies, which also enter 
in the calculation of the wave functions of the outgoing protons and in the 
phase space factors, were taken from [54]. 

3.3 Electromagnetic nucleon current 

The nucleon current 

f N = e* N , f J^ N ,i, (37) 



was modeled by using the ccl current introduced by de Forest, given by the 
operator [55] 

r N {ccl) = {F 1 + kF 2 )Y - {Pl " +P " f) KF 2 , (38) 

2m jv 



with the so-called Dirac and Pauli form factors F 12 (q 2 ) related to the the 
Sachs form factors Ge^Q 2 ), and Gm(q 2 ) [56,57,58], according to 

G E (q 2 ) = F 1 (q 2 ) + j^F2(q 2 ) = F 1 (q 2 ) - nrF 2 (q 2 ), (39) 

G M (q 2 ) = F 1 (q 2 ) + KF 2 (q 2 ), (40) 



where q 2 = (p*j — pf) 2 , and k is the anomalous magnetic moment in units of 
nuclear magnetons (k p = fi p — 1 = 1.792847 for the proton, and K n — /i n — 
— 1.913043 for the neutron). We emphasize that pf, and correspondingly 
q 2 are differential operators acting on the initial and final state nucleon wave 
functions ^nj and ^n,i, and should not be confused naively with correspond- 
ing asymptotic C-number momenta of the nucleon. However, for the sake of 
convenience, we maintain here this notational ambiguity, which is irrelevant 
in the case of free particles described by plane waves. We note that we used 
a completely analogous expansion of the g 2 -dependence of the form factor op- 
erators as it was used by Knoll for the photon propagator in eqns. (8,9,10) in 
our numerical calculations. 

We adopted the dipole formula for the nucleon form factors according to 
GUQ 2 = -<f) = G D (Q 2 ) = (l+ 0J1 Q ey3 ) \ (41) 
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with a magnetic proton form factor given by G P M (Q 2 ) = (1 + k p )G p e (Q 2 ), or 
equivalently 

Ff=^f(P El F p = -^-G p E . (42) 



In the case of the neutron we used the Galster parametrization [59] for the 
(small) electric form factor 

G'k(Q 2 ) = -™„G D (Q 2 )UQ 2 ), (43) 

^ = { 1 + oS^T' (44) 



consistent with the world data [60,61], and G1 I (Q 2 ) = K n G D (Q 2 ), such that 



fin I ~-fin fin fin 

G E + TG M _ G M G E 

1 + r ' 2 K n (l + r) 



An alternative to the ccl current is the cc2 current 

r N {cc2) = F 1 -f + ^ a ^q v F 2 , (46) 
irriN 



which was also used in order to check how the ratio of the cross sections 
behaves for different models of the nuclear current. None of the expressions 
for the current (ccl,cc2) is fully satisfactory and both expressions fail to fulfill 
current conservation, but given the fact that we focus mainly on the electronic 
part of the problem, the simple choices given above provide a satisfactory de- 
scription of the proton current. An advantage of the cc2 choice is the fact that 
integration by parts allows one in a simple way to get rid of the momentum 
operators acting on the nucleon wave functions. 

The asymptotic momenta of the nucleons were calculated from energy and 
momentum conservation, i.e. the energy of the knocked-out nucleons was re- 
duced by their initial binding energies. We also performed calculations with 
enhanced binding energies in order to take into account that the average nu- 
cleon removal energy is higher than the average binding energy. This strategy 
leads to a shift of the quasi-elastic peak to higher energy transfer, but has no 
significant impact on the general conclusions drawn in this paper concerning 
the applicability of the EMA. We found in general that different choices for 
the nucleon current do not affect significantly the relative behavior of PWBA, 
EMA and DWBA calculations. This is not the case for the absolute cross 
section themselves. 
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The accuracy concerning the calculation of cross sections is limited to about 
1% due to the truncation of the Knoll expansion eq. (10) and the finite res- 
olution of the grid that has been used for the modeling of the nucleus. The 
numerical evaluation of transition amplitudes was performed by putting the 
nucleus on a three dimensional cubic grid with a side length of 30 fm and a 
grid spacing of (30/n gr id) fm, and convergence was checked by using different 
side lengths and grid resolutions. The number of necessary grid points ~ n z grid 
is mainly dictated by the wave length of the oscillatory behavior of the matrix 
element eq. (9) ~ e l ( k i- k f-Pf) anc i was j n the range of 70 3 to 130 3 in order to 
ensure an accuracy better than lO^ 1 percent for the values of the integrals. The 
accuracy of the solid angle integration of the (e, e'N) cross section was bet- 
ter than 0.05%. The truncation of the partial wave expansion of the electron 
wave functions [37] was performed such that all partial waves corresponding 
to angular momenta up to \k\ = 70 were taken into account, which guaranteed 
an accurate evaluation of the electron wave functions and corresponding first 
and second order derivatives appearing in the Knoll expansion in the relevant 
nuclear vicinity. 



4 Results 



We first discuss some phenomenological properties of the quasi-elastic peak, 
calculated within the electron PWBA, EMA, and DWBA SPS framework, 
and from an experimental point of view. Within a simplified picture like the 
Fermi gas model, the position of the quasi-elastic peak for (e, e') scattering 
with initial electron energy €i and scattering angle e is expected at 

_Qleak _ C?(l-COs8 e ) 

peak ~ 2m N + €r - ~ m N + (1 - cos0 e ) ei + 1 ] 



where e rem can be interpreted as a phenomenological average removal energy of 
approximately 44 MeV [19]. These 44 MeV include a Coulomb shift of approx- 
imately 8 MeV, such that the physically relevant Coulomb corrected removal 
energy is rather given by ~ 36 MeV, if the comparably small contribution of 
neutrons to the (e, e') cross section is not taken into account. The naive ex- 
pression eq. (47) predicts a difference of the position of the PWBA and EMA 
quasi-elastic peaks of 7.1 MeV for e< = 485 MeV and e = 60°, and 11.7 MeV 
for €i = 310 MeV and 6 e = 143°. A fit of the experimentally measured peak 
positions [17,36] leads to 

ex P - cos g) +m N E 

Upeak rh N + e' i (l-cos#) 1 } 
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with mjv and E as fitting parameters, and e\ = q — V is the effective initial 
electron energy with V = — 20MeV in the present fit. Eq. (48) incorporates 
also the background in the experimental data due to physical mechanisms 
like correlation effects and meson exchange currents which are not included in 
the SPS model. The experimental data are reproduced in a very satisfactory 
manner by eq. (48) for 6 e = 60° with fh N = 721 MeV and E = 12.8 MeV in 
the energy range q = 262. ..645 MeV, and for 6 e = 143° with fn^ = 662 MeV 
and E = 10.0 MeV in the energy range = 140. ..420 MeV. From eq. (48), 
one obtains a difference of the position of the PWBA and EMA quasi-elastic 
peaks of 8.3 MeV for = 485 MeV and 9 e = 60°, and 13.4 MeV for e, t = 310 
MeV and e = 143°. 



The corresponding calculations displayed in Figs. 4 and 5 (Figs. 6 and 7) 
based on the nuclear eikonal model (plane wave approximation) confirm the 
behavior of the quasi-elastic peak position described above in a satisfactory 
way. E.g., the peak shifts are 10.9 MeV (12.8 MeV) for the (310 MeV, 143°) 
kinematics and 9.2 MeV (7.8 MeV) for the (485 MeV, 60°) kinematics. 



A further comment concerns the relative amplitude of the PWBA and EMA 
quasi-elastic peak. An analysis of the Saclay data shows that the maximal 
(e, e') peak cross section a peak = -£^\u, peak scales like a peak ~ e^ 2 ' 9 in the 
initial electron energy region ~ 485 MeV for O e = 60°. Since the EMA cross 
section can be obtained by multiplying the PWBA cross section, obtained 
from effective kinematic quantities, with a focusing factor [(e^ — V)/ei} 2 , a 
ratio of cr™ k BA /a™ k A ~ (504/485) - 9 = 1.035 is expected for the PWBA and 
EMA peak values, if the theoretical model is in acceptable accordance with 
the properties of the physical nucleus. For 6 e = 143°, the scaling behavior is 
~ e~ 2 ' 45 in the energy region ti ~ 310 MeV, predicting an approximate ratio 
v™k BA /°peak A - (329/310) - 45 = 1.027, and ~ e " L9 in the energy region Ci ~ 
224 MeV, such that the peak values of the PWBA and EMA cross section are 
nearly identical, as observed in Fig. 8. We note that the problematic 'staircase- 
like' behavior of the Saclay data, discussed in [21], has only a minor influence 
on the valuations presented above. 



The relative amplitudes and the peak shifts presented in [30] do not exhibit 
the behavior expected from the discussion above in a distinct manner, even 
if one takes into account that a central potential value V(0) ~ —25 MeV 
has been used in [30] for the EMA, whereas an average potential value V = 
—19 MeV has been used for our calculations presented in Figs. 4-8. Especially, 
a pronounced enhancement of the DWBA cross sections with respect to the 
EMA cross sections in [30] is absent in the present calculations, except for the 
case where the final state electron energy is specifically low, as shown in Fig. 
8 below. 
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Fig. 4. The differential cross section for 208 Pb(e,e') scattering for e» = 485 MeV and 
e = 60°, based on the eikonal approximation for the nuclear current. The dotted 
line is the electron PWBA result, the dash-dotted line shows the EMA result with 
an effective potential value V = — 19 MeV, and the solid line displays the DWBA 
calculation. 

For the (485 MeV, 60°) kinematics, a slight shift to the left of the DWBA curve 
with respect to the EMA curve is observed. Still, an EMA-type behavior of 
the DWBA result can be noticed. For lower energy transfer, the influence of 
the optical potential on the nucleons with relatively low energy is stronger and 
becomes less important at higher energy transfer, where the EMA and DWBA 
match almost perfectly. The EMA/DWBA mismatch is virtually absent in the 
case where plane waves are used for the outgoing nucleon wave functions (see 
Figs. 6 and 7), i.e. when the distortion of the final state nucleons is not taken 
into account. From a practical point of view, one may observe that the EMA 
and the DWBA curve match almost perfectly in the energy range 80 MeV < 
6i < 180 MeV, if an effective potential value of —16 MeV is used and the EMA 
results are additionally amplified by 2%. However, from the physical point of 
view, one should argue that an effective potential value of V = —19 MeV is an 
optimal choice, if the energy transfer and the corresponding energy of the final 
state nucleons is high enough such that the impact of the final state interaction 
on the cross described by the optical potential becomes less important. The 
DWBA cross sections are still larger by 1 — 2% than the EMA cross sections in 
this kinematical region, i.e., a slight 'overfocusing' can be observed. Still, one 
should keep in mind that the accuracy of the present calculations is limited 
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Fig. 5. The results corresponding to Fig. 4 for the (310 MeV, e = 143°) kinematics. 

to a similar order of magnitude by the use of the expansion eq. (10). For 
the (310 MeV, 143°) kinematics, the situation is similar, although the impact 
of the final state interaction described by the optical potential in the low 
energy transfer region to the left of the quasi-elastic peak is less pronounced. 
This is also due to the fact that a marginally larger effective potential of 
20 MeV would be appropriate for this kinematics, leading to a better fit of the 
EMA and DWBA for higher energy transfer. Furthermore, the quasi-elastic 
peak is also shifted to higher energy transfer compared to the (485 MeV, 60°) 
kinematics, additionally reducing the influence of the optical potential. 

As expected, in the case of the (224 MeV, 143°) kinematics displayed in Fig. 
8, the EMA starts to fail. The reason is twofold. Firstly, the energy of the 
final state electrons is relatively small, and the focusing effect becomes more 
important [51], such that the DWBA cross section becomes significantly larger 
than the EMA prediction. Secondly, the quasi-elastic peak is shifted towards 
smaller energy transfer, where the optical potential increasingly distorts the 
wave functions of the final state nucleons. The momentum transfer squared is 
still relatively large, and is therefore not a contributing factor to the failure of 
the EMA. Note also that the shifts of the EMA and DWBA peaks with respect 
to the PWBA peak are still of comparable size. Again, in order to describe 
the DWBA result within an EMA framework, one could increase the effective 
potential to a higher value such that the EMA and DWBA curve match better 
in the region of large energy transfer. An other strategy would be to derive 
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Fig. 6. The differential cross section for 208 Pb(e,e') scattering at two different elec- 
tron energies and scattering angles (see also Fig. 6) , calculated by using plane wave 
functions for the final state nucleons (dotted line: PWBA, dash-dotted line: EMA 
with effective potential V = — 19 MeV, solid line: DWBA). The peaks are slightly 
narrower than in the nuclear eikonal model, since effects of the final state interac- 
tions described by the optical potential are absent. 

an optimized effective potential which shifts the maximum of the EMA and 
DWBA peak to the same position and to renormalize the EMA curve by an 
appropriate factor subsequently. 

As a general remark, one may observe that the Coulomb corrections intro- 
duced within the framework of the asymptotic expansion eq. (10) according 
to Knoll depend on derivatives of the electron charge and current densities, 
or, by partial integration, on the nuclear four-current. Since the correspond- 
ing nuclear response is rather smooth (differently from elastic and transition 
form factors to discrete levels), and only leading terms in the expansion are 
relevant, one could expect that the dependence of Coulomb effects on different 
models of the nuclear current is rather weak. 

We finally comment on the electron eikonal calculations (EDWBA) presented 
in [34], which originally seemed to be compatible with the Ohio group results 
in [30]. The PWBA and EDWBA curves in Fig. 3 of [34] show a indeed similar 
behavior as the PWBA and DWBA results in [30]. As in [30], the Coulomb 
corrected cross section is larger at the peak than in the PWBA case, seemingly 
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Fig. 7. The results corresponding to Fig. 6 for the (310 MeV, e = 143°) kinematics. 



in contradiction with the EMA. A typical example from [34] is displayed in 
Fig. 9, where also experimental data taken at Saclay have been included [17]. 
However, there is a simple explanation for this discrepancy. For the EDWBA 
calculations in [34], a focusing value for the electron cross section was used 
which corresponds basically to the central value of the electrostatic potential, 
i.e. V(0) = 25 MeV, which leads to an overestimation of the Coulomb cor- 
rected cross sections. Note that the average momentum of the electrons is well 
described by the eikonal integral and also corresponds also to an average value 
of approximately 19 MeV. 

Therefore, the results presented in [34] can be be corrected by implementing 
the focusing of the electron wave functions obtained from the exact electron 
wave functions in the program originally used in [34], with the result illustrated 
below for Fig. 4 in [34], which is displayed again as Fig. 10 in this paper. 
The figure shows the ratio of the cross sections, calculated in PWBA, with 
the Coulomb corrected cross sections according to the EDWBA and EMA. 
Firstly, the EMA curve has to be adapted to an effective potential of 19 MeV 
instead of 25 MeV. This slightly reduces the ratio (Jpwba/ &ema and moves the 
corresponding dotted curve closer to one (the horizontal line). Secondly, the 
focusing factors of the EDWBA calculations must be corrected. E.g., for = 
485 MeV and e f — 385 MeV, the total focusing factor entering the EDWBA 
for the initial and final state electron originally used in [34] corresponded 
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Fig. 8. PWBA, EMA and DWBA cross section for the (224 MeV, 143°) kinematics 
(nucleon eikonal model, V = -19 MeV). The EMA fails to describe the DWBA 
cross sections reliably due to the low final state electron energy and the influence 
of the optical potential on the final state nucleons for low energy transfer. 
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Fig. 9. Cross sections for (e, e') scattering with = 310 MeV and $ = 143° obtained 
in [34]. The plots displays a similar behavior as the results presented in [30], however, 
the Coulomb corrected eikonal cross sections are too large due an overestimation 
of the electron wave function focusing. A slightly different nuclear model was used 
than in the present work, and the reduced strength due to correlated nucleons had 
been neglected. 



approximately to the numerical value 

(485 + 25) 2 x (385 + 25) 2 



485 2 x 385 2 



1.254. 



(49) 
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Fig. 10. Comparison of Coulomb corrections (e^ = 485 MeV, $ = 60°) for different 
approaches in [34]. For the EMA, an effective potential value of 25 MeV was used. 
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Fig. 11. The same as in Fig. 10, but with corrected focusing in the eikonal calculation 
and an EMA curve obtained by using an average potential V = 18.7 MeV. 

The correct total focusing factor should rather be 
(485 + 19) 2 x (385 + 19) 2 



485 2 x 385 2 



~ 1.189. 



(50) 



Accordingly, the EDWBA cross section has to be reduced by 5.2% at oo = 
100 MeV and the corresponding solid curve moves upwards in the plot. The 
consequences are shown in Fig. 11. Note that for Fig. 11, also locally varying 
focusing factors obtained from exact solutions of the Dirac equation were used 
in conjunction with the eikonal approximation for the phase of the electron 
wave functions. An attempt to calculate corrections to the focusing at some 
distance to the nuclear center has already been presented in [62], which, how- 
ever, does not lead to reliable predictions in the important surface region of the 
nucleus and resulted in the overestimation of the focusing in [34]. For the EMA 
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calculation presented in Fig. 11, a slightly smaller value V = 18.7 MeV than 
in the calculational example above was used, since this value has been deter- 
mined experimentally to be V = 18. 7± 1.5 MeV for 208 Pb [18] by a comparison 
of inclusive electron and positron scattering data. This independent experi- 
mental finding is remarkable, since it also supports that an effective potential 
which is slightly smaller than the naive mean potential V ~ AV(0)/5 ~ 20 
MeV is adequate. One might speculate that the contribution of nucleons near 
the surface of the nucleus to inclusive scattering is enhanced, leading to the 
observed reduction of the effective potential. 



5 Conclusions 

As a basic result of this work we conclude that the EMA with an effective po- 
tential V = —19 MeV is a valid approximation for the description of Coulomb 
distortions in the kinematic region where the momentum transfer squared is 
larger than (300 MeV) 2 (such that the length scale of the exchanged photon is 
smaller than the typical size of a nucleus), the energy of the scattered electron 
is larger than 150 MeV (such that the semiclassical description of the elec- 
tron wave functions in the nuclear vicinity is valid) and the energy transfer 
uj = €i — €f is larger than ~ 140 MeV (such that the distortion of the final state 
nucleon wave functions is moderate). In the case mentioned above, and for a 
typical energy range where the initial electron energy is of the order of some 
hundreds of MeV, the DWBA cross sections are generally larger by 1 — 2% than 
the EMA cross sections. If the energy transfer to is smaller than 140 MeV, the 
DWBA cross section can still be approximated by an EMA calculation with 
a phenomenological effective potential V and a minor amplitude correction. 
However, in such a case one has to rely on the model used to describe the 
nuclear current. 

Other kinematic situations than those presented in this paper have been in- 
vestigated, which showed that the behavior of the forward and backward scat- 
tering angle kinematics presented in this paper is typical. 

The coincidence of EMA and DWBA cross sections is rather impressive. As 
a test, the PWBA and EMA cross section calculations which are based on 
plane Dirac waves for electrons were evaluated by both using Dirac plane 
waves and by performing a limit Z — > in the full Coulomb solver for the 
electron wave functions. Performing such a limit in the DWBA solver leads 
to the same cross sections as obtained from the plane wave calculations, such 
that the observed coincidence served as a numerical test for the two unrelated 
calculational techniques. 

It is instructive to compare Fig. 4 in this work to Fig. 16 in the classic work 
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by Rosenfelder on quasi-elastic electron scattering on nuclei [63]. It should 
be noted that Rosenfelder explicitly mentioned already in his work that the 
effective potential V is close to the mean value of the electrostatic potential 
of the nucleus. However, he then wrote down the explicit expression for the 
central value 



of the potential of a homogeneously charged sphere with radius R, which is 
related to the corresponding mean value by V ~ 4V(0)/5. This minor, but 
not irrelevant misapprehension for the case of the EMA, has propagated in the 
literature since then (see, e.g., [6]). Apart from the fact that the initial electron 
energy is marginally larger and the effective potential used by Rosenfelder is 
given by V = —24.9 MeV, Fig. 16 in [63] can directly be related to the result 
presented in Fig. 4. 

Despite many theoretical and experimental efforts to clarify to what extent a 
suppression of the Coulomb sum rule exists, clear answers to the problem are 
still missing. The origin of the discrepancy between the results presented in 
this paper and the findings of the Ohio group remains hitherto unclear, but 
one may speculate that it could be related to the relativistic nucleon poten- 
tial model used by the Ohio group, which may influence the resulting nuclear 
current (spinor distortion or negative energy contributions) in an unexpected 
way. An overview on different theoretical attempts to implement Coulomb 
corrections in the analysis of experimental data can be found in [64] . Accord- 
ing to a reanalysis of experimental data based on the validity of the EMA, 
which is supported by the present work, Meziani and Morgenstern claimed 
that a suppression of the longitudinal structure function of about 40% exists 
at the effective momentum transfer of 500 MeV/c, and tried to explain the 
suppression by a change of the nucleon properties inside the nuclear medium 
[22]. However, in view of the fact that theoretical questions about the inter- 
pretation of relativistic nucleon potential models persist and further related 
effects like meson exchange currents and correlations have not been taken into 
account in a satisfactory manner in theoretical Coulomb distortion calcula- 
tions up to the present, and since the quality of experimental data may be 
questionable in some cases (see, e.g., [21]), it is advisable to await for the ex- 
pected experimental TJNAF results at the high momentum transfer region, 
in which the relevant correlations become small, whereas Coulomb corrections 
become crucial to extract information about the structure functions and to 
understand the Coulomb sum rule. 
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Appendix A 



The nuclear charge distribution and nucleon optical potentials of heavy nuclei 
are often approximated by the help of a Woods-Saxon distribution. Normal- 
izations and moments of Woods-Saxon distributions are rarely found in the 
literature, however, they can be expressed exactly in terms of polylogarithms. 

The volume integral over the Woods-Saxon distribution with a range f and 
diffusivity a 

p( r ) = (52) 

rv ' l _|_ e (r~r)/a v ' 



can be written with help of the trilogarithm 

oo 

4vr J r 2 p{r)dr = -8na 3 p U 3 {-e f/a ) -> -7rf 3 po {a -> 0), (53) 



where the polylogarithms Li n are defined for \z\ < 1 by 

oo z k 



k=l K 



For \z\ > 1, the analytic continuation of the dilogarithm Li 2 (z) and the trilog- 
arithm Li 3 (z) can be obtained via the functional equations 

Li 2 (z) = -U 2 (l/z) - \ log(-^) 2 - ^, (55) 
Li 3 (z) = Us(l/z) - I \og(-zf - ivr 2 \og(-z). (56) 



For a nucleus with charge number Z and total charge eZ one has 

Po = "dpi' (57) 
The potential energy of an electron in the center of the nucleus is given by 

oo 

= — e J p{r)rdr = —ea 2 p Li 2 (—e r ^ a ) 





aZU 2 (-e f / a ) 
2a Li3(-e*V a )' 



(58) 
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and the rms radius of the distribution is given by 



.2v_io„ 2 Li6(-e f/ ^ 



(r 2 > = 12a 2 , — . (59) 
x 7 Li 3 (-e r / a ) v ' 



This follows also from the general expression for the moments of the Woods- 
Saxon distribution 

oo 

I C-^-dr = -a n+l U n+1 (-e f/a ). (60) 

o 



Note that the corresponding expressions for a homogeneously charged sphere 
can be obtained from the limit a — > and 

re™ 

Li^-e*) ~ for x > n 2 . (61) 

?7>. 



Choosing the typical parameters f = 6.6 fm and diffusivity a = 0.545 fm for the 
electric charge distribution of a 208 Pb nucleus with mass number A = 208 and 
charge number Z = 82, one finds that these values are compatible with an rms 
charge radius of 5.50 fm and a central Coulomb potential of Vq = —25.71 MeV. 
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